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We study optimal control strategies to optimize the relaxation rate towards the fixed point of a 
quantum system in the presence of a non-Markovian dissipative bath. Contrary to naive expecta¬ 
tions that suggest that memory effects might be exploited to improve optimal control effectiveness, 
non-Markovian effects influence the optimal strategy in a non trivial way: we present a necessary 
condition to be satisfied so that the effectiveness of optimal control is enhanced by non-Markovianity 
subject to suitable unitary controls. For illustration, we specialize our findings for the case of the 
dynamics of single qubit amplitude damping channels. The optimal control strategy presented 
here can be used to implement optimal cooling processes in quantum technologies and may have 
implications in quantum thermodynamics when assessing the efficiency of thermal micro-machines. 


I. INTRODUCTION 

Controlling quantum systems by using time-dependent 
fields [lj is of primary importance in different branches of 
science, ranging from chemical reactions @,|,NMR a, 
molecular physics 0 to the emergent quantum technolo¬ 
gies @-[ 8 ]. Investigations on optimal control of open 
quantum systems mostly focus on memoryless environ¬ 
ments [9Mll1| and specifically on those situations where 
the reduced dynamics can be described by a Markovian 
master equation of the Lindblad form [12|. In this con¬ 
text, optimal control applications to open quantum sys¬ 
tems have been explored in different settings iTcj - jlsIl 
and recently the ultimate limits to optimal control dic¬ 
tated bv quantum mechanics in closed and open sys¬ 
tems jl9l - l22l ] and the complexity of dealing with many- 
body systems flPl Idbl have been determined. Time- 
optimal quantum control has been extensively discussed 
for one qubit systems in a dissipative environment (9jj26| 
and the optimal relaxation times determined in 1271 ]. 
These studies might have both fundamental and prac¬ 
tical applications, for example in assessin g th e ultimate 
efficiency of quantum thermal machines [281 ] . or to im¬ 
plement fast cooling schemes which have already proven 
to be advantageous HHH. 

However, introducing a Markovian approximation re¬ 
quires some constraints on system and environment, 
which may not be valid in general [3lll32| . Consequently, 
incorporating non-Markovian (NM) effects of the envi¬ 
ronment, in a sense that will be defined more precisely 
below, might be a necessity in a many experimental situ¬ 
ations. Recently, the possible influence of memory effects 
on the ortho gona lity catastrophe [33], on quantum speed 
of evolution [TldT and on quantum control [3^, [36| have 
been analyzed. Here, we present a study of the optimal 
control strategies to manipulate quantum systems in the 
presence of NM dissipative baths and compare the per¬ 
formance of optimal control with the case of operating 
subject to a Markovian (M) environment. 

Intuitively, the absence of memory effects in the dy¬ 


namics of open quantum systems is linked to the possi¬ 
bility of identifying well separated time scales in the evo¬ 
lution of system and environment. Recently, a number of 
proposals have been put forward to quantitatively char¬ 
acterize this effect in terms of explicit non-Markovianity 
measures [37l440l | . In this light, one can define an evo¬ 
lution to be Markovian if described by a quantum dy¬ 
namical semigroup (time-homogeneous Lindbladian evo¬ 
lution) [II} , which would be the traditional Markovianity 
considered in most previous work on open system control. 
However, other definitions encompass this as a special 
case while allowing for more general, non-homogeneous 
generators, albeit still ensuring the divisibility of the as¬ 
sociated dynamical map and the unidirectionality of the 
system-environment information flow, and therefore the 
absence of memory effects in the dynamics of the sys¬ 
tem [42U451 ]. Relevant for our analysis is the definition 
of Markovian evolution in terms of the divisibility of the 
associated dynamical map [H]. When the dynamics is 
parametrized using a time-local master equation, the re¬ 
quirement of trace and hermiticity preservation, yields a 
generator of the form 

Ps(t) = -i[H s (t),p(t)]+C(t)(p(t)) 

= -i[ H s(t),p(t)] + ^lk(t)(A k (t)p(t)Al(t) ( 1 ) 
k 

- ^{A{(t)A k (t)p(t)}), 

where C(t ) is a time dependent Lindblad superopera¬ 
tor, 7 k {t) are generalized (i.e. not necessarily positive) 
decay rates, the 7fi(f)’s form an orthonormal basis for 
the operators for the system, see e.g. Ref. [47| (here¬ 
after h has been set equal to one for convenience) and 
H s (t ) is the effective Hamiltonian acting on the system. 
Equation © generalizes the familiar Lindbladian struc¬ 
ture to include NM effects while maintaining a time-local 
structure. However, apart from same special cases, it 
not known which are the conditions which H s (t), A k (t ), 
and 7 *, ( t) have to satisfy in order to guarantee Complete 
Positivity (CP) [47l45l1 | . i.e. the fundamental prerequi- 
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site which under fairly general assumptions is needed 
to describe a proper quantum evolution [3ll . [32| . In 
what follows we will focus on a simplified scenario where 
the 7 fe(t)’s either are null or coincide with an assigned 
function 7 (f), and where the Ak(t )'s are explicitly time- 
independent. Accordingly, in the absence of any control 
Hamiltonian applied during the course of the evolution, 
we assume a dynamical evolution described by the equa¬ 
tion 


p(t) = 7 (t)£(p(f)) , 


( 2 ) 


where £ is a (time-independent) Lindblad generator char¬ 
acterized by having a unique fixed point pf p (i.e. C(p) = 
0 iff p = Pf p ). For this model, in the absence of any 
Hamiltonian term (i.e. H s (t) = 0) CP over a time inter¬ 
val [0, T] is guaranteed when [48| 

> 0 , Vf € [0,T] , (3) 

while divisibility (i.e, Markovianity) is tantamount to the 
positivity of the single decay rate at all times [13 : if 
there exists a time interval where 7 (f) becomes negative, 
the ensuing dynamics is no longer divisible and the evo¬ 
lution is NM. In this context we will assume a control 
Hamiltonian H s (t ) to represent time-localized infinitely 
strong pulses, which induce instantaneous unitary trans¬ 
formations at specific control times. This corresponds to 
writing H s (t) = JT <5(t — where 0^ are time 

independent operators which act impulsively on the sys¬ 
tem at t = tj ( S(t ) being the Dirac delta-function), at 
which instants one can neglect the contribution from the 
non-unitary part, and represent the master equation by 
p{t) ~ — i [H s (t), p{t)}. Therefore the resulting dynamics 
is described by a sequence of free evolutions induced by 
the noise over the intervals t £ [tj, tj + 1 ] interweaved with 

unitary rotations Uj = exp[—* 0 ^], i.e. 

p{t) = Ufin O Vj O Uj O Dj- 1 o • • • 

■■•oW 1 o’D 0 oW m (p(0)), (4) 



where = exp f’ i+1 C(t) , U{- ■ ■) = [/(• • -)W and “o” 
is the composition of super-operators. When only two 
control pulses are applied (the first Ui n at the very begin¬ 
ning and the second Ufi n at the very end of the temporal 
evolution), the non-unitary evolution is described by Eq. 
© and CP of the trajectory © is automatically guaran¬ 
teed by Eq. ©, the scenario corresponding to the real¬ 
istic case where one acts on the system with very strong 
control pulses at the state preparation stage and imme¬ 
diately before the measuring stage. When more Uj 's are 
present, the situation however becomes more complex. 
There is no clear physical prescription which one can fol¬ 
low to impose the associated dynamics on the system at 
least when the dissipative evolution is assumed to be NM. 
Consider for instance the case Ufi n o'DioUioT>ooUi rl (p(0)) 
where U\ is a non trivial unitary. Even admitting that 
the latter is enforced by applying at time t\ a strong 


instantaneous control pulse, there is absolutely no clear 
evidence that the open dynamics for t>t\ should be still 
described by the same generalized Lindbladian 7 {t)C, the 
system environment being highly sensitive to whatever 
the system itself has experienced in its previous history. 

We note that in a more realistic scenario, any control 
pulse will have a non-zero width St in time. Clearly, a 
sufficiently large St can invalidate the assumption of ap¬ 
plying control pulses only at the very beginning and the 
very end, thereby modifying the dynamics significantly as 
described above. However, one can expect the assump¬ 
tion of instantaneous pulses to be valid as long as St is 
negligible compared to the time scale associated with the 
dynamics in absence of any control. 

Keeping in mind the above limitations, in this work we 
attempt for the first time a systematic study of optimal 
control protocols which would allow one to speed up the 
driving of a generic (but known) initial state p( 0 ) toward 
the fixed point pj p of the bare dissipative evolution for 
the model of Eq. m which explicitly includes NM effects. 
We arrive at the quantum speed limit times when appli¬ 
cation of only two control pulses Ui n and Ufi n , at initial 
and final times respectively, is enough to follow the op¬ 
timal trajectory. On the other hand, we present lower 
bounds to the same when optimal control strategy de¬ 
mands unitary pulses at intermediate times as well. We 
show that the efficiency of optimal control protocols is 
not determined by the M/NM divide alone but it depends 
drastically on the behaviour of the NM channel: if the 
system displays NM behavior before reaching the fixed 
point for the first time, NM effects might be exploited 
to obtain an increased optimal control efficiency as com¬ 
pared to the M scenario. On the contrary, NM effects 
are detrimental to the optimal control effectiveness if in¬ 
formation back-flow occurs only after the system reaches 
the fixed point (see Fig. [TJ. These results are valid ir¬ 
respective of the detailed description of the system, i.e. 
its dimension, Hamiltonian, control field, or the explicit 
form of the dissipative bath. 


II. THE MODEL 

The divisibility measure for the model Eq. m is equiv¬ 
alent to the characterization of memory effects by means 
of the time evolution of the trace distance [ 53 . This 
provides an intuitive characterization of the presence of 
memory effects in terms of a temporary increase in the 
distinguishability of quantum states as a result of an in¬ 
formation back-flow from the system and into the envi¬ 
ronment that is absent when the evolution is divisible 
[53| . As a result, a divisible evolution for which the sin¬ 
gle decay rate 7 (t) > 0 at all times will exhibit a mono¬ 
tonic decrease of the trace distance of any input state 
towards a (assumed to be unique) fixed point pf p of the 
Lindblad generator C [54|. On the contrary, as illus¬ 
trated in Fig. [I] the behaviour of the trace distance can 
be non-monotonic when the dynamics is NM. In this case, 
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FIG. 1: (Color online) Schematic diagram of NM dynamics in 
Class A (red line) and Class B (black line). The instantaneous 
trace distance d(t) = || p(t) — pf p || starts increasing only after 
the system reaches the fixed point when | 7 t| —» oo in case 
of Class A, while it shows oscillatory behavior even before 
it reaches the fixed point in case of Class B. In comparison 
dynamics for a M channel is shown by the blue line where d(t) 
decreases monotonically and assymptotically to d(t — ¥ oo) = 
0. The speedup obtained by M dynamics is always bigger than 
that obtained for the NM one, i.e., Rm/Rnm > 1 in case the 
NM evolution is of Class A, while the Markovian limit can be 
surpassed by NM of class B. 


there exist time interval(s) where 7 (t) becomes negative. 
Denoting by d(t) = || p(t) — pj p \\ the trace distance be¬ 
tween p(t) and the fixed point pf p , it straightforwardly 
follows that d(t) < 0 Vi in the M limit. Looking at this 
quantity one can classify NM dynamics into two distinct 
classes (see Fig. [TJ: the first one (Class A) is defined by 
those dynamics where the system reaches the fixed point 
at time Tp before q(t) changes sign i.e., 7 (f) > 0 and 
d(t) < 0 for 0 < t < Tp. In this case, the NM dynamics 
reaches the fixed point pf p and then start to oscillate. 
On the other hand, Class B dynamics is characterized 
by 7 (f) that changes sign (and correspondingly d(t) > 0) 
at some time t < Tp, that is the solutions of the equa¬ 
tion 7 (t s ) = 0 are such that t, s < Tp for at least one 
s. In contrast, in the M dynamics d(t) always decreases 
monotonically and asymptotically to d{t —> 00 ) = 0 . 

NM channels of class A/B arise from different physical 
implementations. As an illustration, the damped Jaynes 
Cummings model exemplifies a Class A dynamics. Here 
a qubit is coupled to a single cavity mode which in turn is 
coupled to a reservoir consisting of harmonic oscillators in 
the vacuum state (see Eq. © ) 0 , fPCl ICTr . On the other 
hand, dynamics similar to Class B can arise for example 
in a two level system in contact with an environment 
made of another two level system, as realized recently in 
an experimental demonstration of NM dynamics [HU. 

As we will see hereafter, the difference between Class 
A and B appears to drastically affect the performance 
of any possible optimal control strategy to improve the 
speed of relaxation of the system towards the fixed point. 


We assume full knowledge of the initial state and we 
allow for an error tolerance of 0 < e 1 , considering that 
the target is reached whenever the condition |d(t)| < e 
is satisfied. To obtain a lower bound on the minimum 
time Tqsl needed to fulfill such constraint we restrict 
our analysis to the ideal limit of infinite control which 
allows us to carry out any unitary transformation instan¬ 
taneously along the lines of (and with all the limitations 
associated with) the formalism detailed in Eq. ©. In 
the limit of infinite control an important role is played 
by the Casimir invariants r, (j = 2,3for a N 
level system). The Casimir invariants of a state p are 
related to the trace invariants Tr(p') (j = 2, 3,..., N) and 
they cannot be altered by unitary transformation alone 
fiol |59 ]. For example, a two level system has a single 
Casimir invariant its purity P = Tr (p 2 ), which remains 
unchanged under any unitary transformation. Conse¬ 
quently, any optimal strategy with the controls restricted 
to unitary transformations only, would be to reach a state 
p characterired by all Casimir invariants same as those 
of pf p in the minimum possible time. Following this we 
can apply a unitary pulse to reach the fixed point instan¬ 
taneously. Clearly, any constrained control will at most 
be as efficient as the results we present hereafter, based 
on the analysis we have presented previously for the case 
of M dynamics @, 0 . 

In what follows, we will analyse Class A and B channels 
independently. 

Class A: As shown in Fig. [TJ in the NM regime d(t) 
goes to zero at t = Tp when p = pj p and C{pf p ) = 0. 
At the same time we expect | 7 (f)| —> 00 at t ~ Tp in 
order to have finite p(t) = 'y(t)C(p(t)) even for C(p(t)) ~ 
C(pf p ) = 0 , as is required for a non-monotonic d of the 
form shown in Fig. [TJ Notice that 7 (f) and hence the time 
t = Tp at which 7 (f) —> 00 are in general independent of 
p. Consequently any optimal control protocol which in¬ 
volves unitary transformation of pit) generated by H s (f) 
at earlier times t < Tp followed by non-unitary relaxation 
to pf p is expected to be ineffective in this case and we 
have Tqsl = Tp. That is, the gain (or efficiency) of opti¬ 
mal trajectory in the NM class is R^ M = Tf/Tqsl = 1- 
One can easily see Tp/Tqsl = 1 implies absence of any 
speed up, whereas any advantage one gains by optimal 
control can be quantified by Tp/Tqsl > 1- On the other 
hand in the M limit 7 (f) = 70 is finite and constant, and 
the system relaxes asymptotically to the fixed point in 
the absence of any control. In this case we introduce an 
error tolerance Kl, such that we say the target state 
is reached if \d(Tp)\ < e. Clearly, Tp increases with de¬ 
creasing e diverging to Tp —y 00 in the limit e —► 0, as can 
be expected for finite 70 . Therefore the above argument 
of 17 (f) | —> 00 at t ~Tp does not apply in this case and 
in general one can expect the time of evolution to depend 
on the initial state. Consequently the quantum speed up 
ratio Rfj can exceed R^m ~ 1 ; as is explicitly derived 
below in the case of a two level system in presence of an 
amplitude damping channel. Similar arguments apply 
also in the case when an additional unitary transforma- 
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tion is needed at the end of the evolution to reach pf pi 
where — > oo for e — > 0 [13 • 

Our above result > R^m can be expected 

to be valid in a more generic scenario with p(t) = 
lk(t)Ck(p{t)) as well, where not all y^’s 0) are 
same, Ck s are the time independent Lindblad generators 
and the unique fixed point pj p is defined by Ck(pfp) = 0 
for all k. In this case at least one of the can be ex¬ 
pected to diverge at time t = Tp in order to ensure Class 
A NM dynamics as shown in Fig m, thus making any op¬ 
timal control ineffective as detailed above. We note that 
one can have dynamics with time dependent Lindblad 
generators and uncontrollable drift Hamiltonians acting 
on the system during the course of the evolution, in ad¬ 
dition to the instantaneous control pulses, as well. The 
drift Hamiltonians can be expected to modify the Lind¬ 
blad generators thus making the problem more complex; 
however the analysis in this case is beyond the scope of 
our present work. 

Class B: Here we focus on systems of Class B where 
as already mentioned y(t) changes sign for t s < Tp with 
s = 1 ,N S . Clearly, in this case i'y(t)| does not nec¬ 
essarily diverge for any t. Consequently the arguments 
presented above for class A fails to hold any longer and 
the time of relaxation to the fixed point can in general be 
expected to depend on p(t) (and hence on H s ). Further¬ 
more, it might be possible to exploit the NM effects such 
that even though d(t) > 0 for t\ < t < t 2 one can, by ap¬ 
plication of optimal control, make sure that Tj(t) > 0 and 
maximum V t, j (where we have assumed Tj(t = 0 ) < Fjf 
V j and Fjf denotes the jth Casimir invariants for the 
fixed point Pf p ). This presents the possibility of exploit¬ 
ing NM effects to achieve better control as opposed to 
the M dynamics, as is presented below for the case of a 
two level system in the presence of an amplitude damp¬ 
ing channel. However we stress that this is not a general 
result and explicit examples can be constructed where 
this is actually not true. 


A. Generalized amplitude damping channel 

Let us now analyze in detail the generic formalism out¬ 
lined above for the specific case of a two level system in 
contact with NM amplitude damping channels of the two 
classes introduced before. 

We consider the non-unitary dissipative dynamics de¬ 
scribed by the time local master equation acting on a 2 x 2 
reduced density matrix p(t) of a qubit and we consider 
the time independent Linbladian C given by 

C(p{t)) = Ci{p{t)) +e /3 £ 2 (p(i)), 

A (p(t)) = ^cr + p(t)cr_ - i{cr_cr + ,p(t)}^ , 

£2 (p{t)) = [a-p(t)(T + - i{cr+CT_,p(t)}^ , ( 5 ) 

with a± being the raising/lowering qubit operators and 


1 //? gives the temperature of the bath. The system evo¬ 
lution is given by Eq. [ 2 j and we will focus on two dif¬ 
ferent functional dependence of the parameter 7 (t) cor¬ 
responding to the Class A and B dynamics. We will 
analyze the system evolution following the Bloch vector 
r representing the state p = (J + r.a) /2 inside the Bloch 
sphere, where the unitary part of the dynamics gener¬ 
ated by H s induce rotations, thus preserving the purity 
P = (1 + |r| 2 )/ 2 . In contrast, in general the action of the 
noise is expected to modify the purity as well. 

Class A: An example of this class of dynamics is ob¬ 
tained under the assumption 


7 (t) = 


2A70 sinh 


g cosh % + A sinh ^ 


g = VA 2 - 270A. 


( 6 ) 


In the above expression A and 70 are two positive con¬ 
stants whose ratio determines the bath behavior: A > 270 
corresponds to a M bath, whereas g becomes imaginary 
in the limit 70 > A /2 resulting in NM dynamics. In the 
NM limit of 70 3 > A, 1 the bath time scale is determined 
by the product A70 and is independent of the specific 
form of the super-operator C. It can be easily seen that 
7 (t) increases monotonically from 0 to 7 (t) — > 00 at 


lim Tp « , 

7o/a->oo V2A70 


( 7 ) 


where Tp is independent of the initial state and the sys¬ 
tem reaches the fixed point when y(£) diverges. With this 
choice of y(£) in Eq. © the time scale is given by 
in the NM limit 70 3 > A, while 7 (t) « 70 sets the time 
scale in the M limit A 70. Therefore the time taken 
to reach the fixed point can be expected to decrease as 
1/VWo i n the NM limit while it scales as ~ I/70 in the 
M limit. 

As mentioned above, this form of 7 (t) arises in the 
damped Jaynes-Cummings model at absolute zero tem¬ 
perature, where one considers only a single excitation in 
the qubit-cavity system and Eq. © reduces to C(p(t)) = 
£■2 However, here we consider a phenomenological 

form of the Lindblad generator © with arbitrary /3 to 
show the generality of our results. In this context the pa¬ 
rameter A in 7(£) denotes the spectral width of the cou¬ 
pling to the reservoir, while 70 characterizes the strength 
of the coupling. 

In the absence of any control the qubit relaxes to 
a fixed point pj p characterized by the Bloch vector 

^0,0, ^ and the optimal control we analyze here aims 

to accelerate the relaxation towards this state with un¬ 
constrained unitary control. Following the strategy pro¬ 
posed above, we look for the extremal speed v = dP/dt 
of purity change for every r. For this model, the speed 
of change of purity is given by v(r, 9 ,t ) = —7(£)(e' S — 


1 )r 


cos 9 


2r fi 


(l + cos 2 0) : note that a positive (neg¬ 


ative) v denotes increasing (decreasing) purity. The two 
strategies differ slightly in case of cooling or heating (i.e. 
the final purity is lower or higher than the initial one); 
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FIG. 2: (Color online) (a) Parametric plot showing variation 
of time Tq£ l of reaching the fixed point with A and 70 for 
/3 = 2, r, = 0.5 and e = 0.01. The Markovian (M) and non- 
Markovian (NM) regions are separated by the blue line on the 
A — 70 plane, (b) Plot showing variation of quantum speed up 
ratio R a with A and 70 for (3 = 2, r x % = 0.3, r y i = 0, r z i = 0.4 
and e = 0.01. R A saturates to Rm « 2 (Rnm ~ 1) in the 
extreme M (NM) limit. 


but both cases correspond to applying unitary rotations 
at the beginning and at the end (for heating) of the 
dynamical evolution, thus yielding a trajectory of the 
form Eq. 0} which is fully compatible with the CP re¬ 
quirement and which doesn’t pose any problem in terms 
of physical implementation (see discussion in Sec. HT1) . 
Specifically we need to apply unitary control so that the 
system evolves along 9 = n till the final purity is reached 
in the case of cooling, while 9 = 0 is the optimal path in 
the case of heating [271 ], in agreement with a recent work 
on quantum speed limit in open quantum systems [HJ 
(see Appendix IIVI for details). Our analysis clearly shows 
that TggL decreases with increasing A as ~ l/\/A for 
small A, finally saturating to A independent constant 
values in the M limit (A 1). However, it would be 
misleading to conclude about the role of Markovianity 
on Tqsl from this alone, since both T l Qg l L and Tqsl de¬ 
crease with increasing 70 as well. In particular, they scale 
as ~ I /70 in the M limit of small 70 , while the scaling 
changes to l/^To for A /70 —1 0. Indeed, the behavior 
depends more on the specific path in the ( 70 , A) plane 
rather than whether the system is M or NM (see Fig.©i). 
However, if one analyzes the speedup obtained by means 
of optimal control strategies, the scenario changes: in 
this case (as shown in Fig. [^Jo) the ratio R A = Tp/Tg™ l L 
clearly distinguishes between M and NM dynamics with 
typical limiting values given by lim € _>o,A/ 7 0 ->-oo —1 2 

and lim e _ > o,A/ 7 0 -s-o Rnm 1 ( see Appendix for details). 
Finally, one can show that a control pulse can be con¬ 
sidered to be instantaneous as long as its width in time 
St -C 1/ [ 70 (e^ + 1)] in the M limit of A 70 , while in 
the NM limit A 70 one has St -C l/i/A7o(e^+T)- 
Class B: Finally, we investigate a particular case be¬ 
longing to the class B dynamics and compare it to the 
previous case. As in the previous case we consider a time 



FIG. 3: (Color online) (a) Plot showing the variation of 
the ratio of the gains Rnm/R-m of the optimal trajectory 
as a function of fl in case of Class B (cooling) with 7 (t) = 
exp(— t) cos(flt), P = 2 , e = 0.01 and initial state given by 
r x i = 0.3, r y i = 0, r zi = 0.4. Clearly gain Rm (= 3.4) in 
the M limit £7 = 0, shown by the green dot, is less than 
that (Rnm) in the NM limit £1 > 0. (b) Schematic diagram 
showing optimal path in the x — z plane of the Bloch sphere 
(red curve) in case of cooling a qubit in Class B (Eq. (fSjl). 
when we start from an arbitrary state pi. The fixed point p/ 
is denoted by brown star. 

evolution described by a master equation Eq. © with 
C given by Eq. ©; however for our present purpose we 
formulate a 7 (f) given by 

7 (t) = e~^cos(fit) , ( 8 ) 

with C, n being two positive constants satisfying the CP 
condition Eq. ©. With this choice, in the absence of the 
control Hamiltonian, NM effects manifest themselves for 
(2 n + 1)7 t/2 < f It < (2 n + 3 )tt/2 for integer n > 0 as 
7 (t) changes sign at f It = (2 n + 1)7t/ 2, simultaneously 
altering the sign of d(t) to d(t) > 0. With a proper 
choice of parameters one can make 7 (t) (and hence d(t)) 
exhibit oscillatory dynamics for 1 > d > 0. As for the 
previous example, also in this case the extremals of v 
are independent of 7 (f) and determined by C(p(t)) only. 
Therefore they occur at exactly the same points as for 
class A ©, i.e., at 9 = 0 , 7 r and arccos (r/rf). In this 
case an instantaneous pulse would correspond to its time 
width St < min{l/£, 1 /Cl, 1/ (e^ + l)}. 

The unconstrained optimal strategies are now modi¬ 
fied as follows. In the case of cooling, the optimal strat¬ 
egy is to follow the path 9 = tt for 0 < t < Tqq 1 l , 
during which time the purity increases monotonically, 
where we have assumed the system reaches the target 
at t = Tqsl < 7r/(2fi) for simplicity. Consequently the 
optimal strategy demands a single pulse at time t — 0 
only to make 9 = 77 which corresponds to an evolution 
operator of the form Ufi n o D o Ui n . Clearly, this evo¬ 
lution is CPT as already discussed in section Q] with D 
depending on and C (see Eqs. © and ©) and is 
thus possible to implement physically. In Fig. ©a) we 
report the speedup obtained by such optimal strategy 
for different values of f 1 in Eq. ©, where we have taken 
e = 0.01 large enough so that flTgg’/ < n/2. In this 
case we arrive at the M limit by setting fl = 0; as can 
be clearly seen the speed up in the NM limit is such that 
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R^m > R-m j ^ ^ > 0 showing that there exist scenarios 
where NM effects can be exploited to improve the con¬ 
trol effectiveness. On the contrary, if the system does 
not reach the target for MTqsl < 7 r/ 2 , the optimal strat¬ 
egy changes: at fit = ir/2, 7 (t) and hence v change sign, 
leading to decrease (increase) of purity for 9 = ir (9 = 0). 
Interestingly, we can take advantage of this effect by mak¬ 
ing 6 = 0 at t = 7 r/( 2 f1), where v exhibits a maximum for 
7 (t) < 0. As mentioned before, application of a unitary 
pulse during the course of an evolution may change the 
form of 7 (t) and C (see Eq. <(4])). However, for simplicity 
let us assume v changes sign at t = (2 n + l)7r/2 and as¬ 
sumes extremum values at 9 = 0 , 7 r and arccos (— r/rf p ) 
even in presence of unitary control. One can easily extend 
our analysis to a more generic case where the simplifying 
assumptions do not hold by following the path of max¬ 
imum (minimum) v for cooling (heating). Let us first 
consider the case where the system reaches the target 
r z = rf — e at t = Tqsl < 37 r/( 2 f l). In such a scenario, 
as depicted in Fig. [3] (b), we let the system evolve freely 
for 7t/(2H) < t < Tqs 1 2 3 4 5 6 7 l , following which we take the sys¬ 
tem to 9 = tt and r z = —rg+c, thus obtaining the desired 
goal. Clearly, an optimal path exists in case of the class B 
non-Markovian channel, which if possible to be followed 
by application of suitable unitary controls, helps in cool¬ 
ing and in particular might make it possible to reach the 
fixed point in finite time (if e = 0). Generalization to 
the case QTq° g l L > 3n/2 where multiple 7 r rotations are 
needed is straightforward. However, we emphasize that 
the strategy presented above for > n/2 follows 

an evolution operator of the form Eq. (U) with unitary 
pulses applied at intermediate times. Consequently our 
analysis gives a lower bound to Tqsl only, achieved by 
following the optimal path shown in Fig. [3jb) , for which 
at present we do not have any implementation strategy. 

Finally, we address the problem of heating the system 
in the shortest possible time, which amounts to minimiz¬ 
ing v V t. Therefore, the optimal path dictates to set 
9 = 0 at t = 0 and then let it evolve freely till f It = 7r/2, 
where 7 (t) and hence v change sign. Unfortunately in 


contrast to the cooling problem, now v > 0 \/9 making it 
impossible to take advantage of the NM effects to accel¬ 
erate the evolution. However, even in this case one can 
always minimize the unwanted effect of information back- 
flow for 7 (t) < 0 by increasing 9 to 9 = arccos(— rj p /r) 
where v has a minimum. 


III. CONCLUSION 

We have studied the effectiveness of unconstrained 
optimal control of a generic quantum system in the 
presence of a non-Markovian dissipative bath. Contrary 
to common expectations, the speedup does not crucially 
depend on the Markovian versus non-Markovian divide, 
but rather on the specific details of the non-Markovian 
evolution. We showed that the speed up drastically 
depends on whether the system dynamics is monotonic 
or not before reaching the fixed point for the first 
time, as determined by the trace distance to the fixed 
point (class A and class B dynamics respectively). 
Indeed, in the former case, the speed up obtained 
via optimal control is always higher in presence of 
a Markovian bath as compared to a non-Markovian 
one, while the reverse can be true in the latter case. 
Finally, we have presented some specific examples of 
these findings for the case of a two level system subject 
to an amplitude damping channel. Note that, in the 
more realistic scenario where one can apply control 
pulses of finite strength only, the presented results serve 
as theoretical bounds to the optimal control effectiveness. 
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IV. APPENDIX 


A. Class A optimal strategy 


Analysis of v{r,6,t) 
1 )r cos 0 + (l + c°s 2 0) 


-70) (e^ - 


shows the optimal strat¬ 


egy in case of cooling is to apply a unitary pulse at 
t = 0 so as to rotate r to 9 = n. Following this we 
switch off the control and allow the qubit to relax by 
the application of the dissipative bath alone for a time 
t = Tqq 1 l , till it reaches r(TQg l L ) = r/ p — e. In contrast, 
while considering the problem of minimizing the time 
taken to heat the qubit, the optimal strategy is to first 
rotate the Bloch vector to 8 = 0. As before, we then 
turn off the unitary control and let it relax till it reaches 
( 0 , 0 , rf p ), following which we apply a second unitary 
pulse to take the system to the fixed point r/ p [27J ■ 

We use the optimal strategy formalism presented above 
to arrive at the minimum time Tq needed to cool the 
system in the different limits. The time for cooling in the 
M limit A /70 —> 00 is given by 


lim 

A/ 70 —>-oo 


rpCOOl 

1 QSL 


1 

7o(l + eP) 


In 


r fp - n 
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(9) 






whereas the same in the NM limit is 


lim T, 

A/70—>-0 


cool 

QSL 


2 

A70 


r fp 


[2(exp(/3) + l)] 


( 10 ) 


On the other hand, following the optimal strategy to heat 
the qubit one gets 


where we have assumed r X i ^$> e, which is typically the 
case for e —>■ 0. In the case of cooling, as can be seen 
from Eq. ([9]) and Eq. da, R A t rs 2 in the M limit as 
long as r X i , ( rf p —ri)^>e and e —> 0. On the other hand, 
the NM limit of A /70 —> 0 yields 


lim Tkfl « —- 3 - 

a/7o->oo 7o( 1 + e^ 3 ) 


In 




/ 2 

' n ( e \ 2 /( exp ( /3 )+ 1 )" 

'rfp + n 
. 2 r fp + e. 

T f m 

( 11 ) 

V A 70 

2 V r *o J 


(14) 


in the M limit, while in the NM limit it is 


lim Tpfft 
A/7o->0 QSL 



/ 2 rf p + e\ 2(ex P o)+1) 
\ri+r fp ) 


As can be easily seen, one needs to consider a non-zero e 
in order to keep the time of cooling finite in the M limit 
while we can set it exactly to 0 in the other cases and yet 
reach the fixed point in finite time. 

In contrast to the results derived above, the advantage 
one gains by application of optimal control presents a 
completely different picture. As before, one can under¬ 
stand this from the the quantum speed-up ratio R A = 
Tf/Tqsl- la the M limit A /70 —>■ 00 , 7 (f) « 70 one gets 


T f 


2 

70 (1 + e&) 


In 


6 


(13) 


thus reducing the gain to R^ M ~ 1 in the limit of e —> 0 
(see Fig. (EJr)). In case of heating we have 


lim R 

e ->0 
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Tfp + Ti 

2 r fp 
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in the M limit, while the same in the NM limit is 


lim Rnm 


7t/2 



r in 

cos -1 

/ 2r fp \ 2(exp(/3) + l) 
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(16) 


which again implies the gain is much higher in the M 
limit for e —> 0 . 






























